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ABSTRACT 

In this paper we present a new method to estimate a foreground cleaned Cosmic Microwave Background (CMB) 
map at a resolution of 1 ° by minimizing the non-Gaussian properties of the cleaned map which arise dominantly 
due to diffuse foreground emission components from the Milky Way. We employ simple kurtosis statistic as the 
measure of non-Gaussian properties and perform a linear combination of 5 frequency maps provided by Wilkinson 
Microwave Anisotropy Probe (WMAP) in its 7 year data release in such a way that the cleaned map has a minimum 
kurtosis which leads to a non-Gaussianity minimized, foreground cleaned CMB map. We validate the method 
by performing Monte-Carlo simulations. To minimize any residual foreground contamination from the cleaned 
map we flag out region near the galactic plane based upon results from simulations. Outside the masked region 
our new estimate of CMB map matches well with the WMAP's ILC map. A simple pseudo-C/ based CMB TT 
power spectrum derived from the non-gaussianity minimized map reproduces the earlier results of WMAP's power 
spectrum. An important advantage of the method is that it does not introduce any negative bias in angular power 
spectrum in low multipole regime, unlike usual ILC method. Comparing our results with the previously published 
results we argue that CMB results are robust with respect to specific foreground removal algorithms employed. 

Subject headings: cosmic background radiation — cosmology: observations — diffuse radiation 



1. INTRODUCTION 

The Cosmic Microwave Background (CMB) has conceiv- 
ably become the finest tool so far to probe the physics of the 
early universe. The problem of isolating a clean CMB sig- 
nal from contaminations originating from Milky Way is of pri- 
mary importance to a cosmologist. According to the slowly 
rolling single scalar field inflationary scenario the primordial 
perturbations follow Gaussianity to a very good approxima- 
tion, a ny non Gau s sian p r operties pr e dicted t o be a small 
effect dAllen et al.l (11987): lFalk etal.1 (|1992|): iGangui et all 
( 11994 : lAcquavivaet all (l2003h : iMaldacenal ( I200 3U Since 
CMB anisotropics are directly related to these perturbations 
via the spacetime metric one expec ts CMB to follow G aus- 
sian distribution. It was shown by 'M unshi et al.l d 19951) that 
non-Gaussian effect introduced in CMB after decoupli ng is 
also small. It has been shown by Komatsu et al. ( 2003) that 
CMB data observed by WMAP satellite (Bennett et al. (1997) ) 
is consistent with primordial Gaussian fluctuation. Contrary 
to these, diffuse foreground components, originating from 
our own galaxy exhibit highly non-Gaussian properties due 
to non-linearities involved in the physical processes during 
their origin. Based upon simple assumption of pure Gaus- 
sian properties of CMB and non-Gaussian nature of diffuse 
foregrounds, in this paper, we propose a new foreground re- 
moval method from CMB maps. A detailed account on sev- 
eral other foregro und removal method s may be found in ex- 
isting Uteratures ("Br andt et alj (11994 : iTegmark & Efstathioul 

(2003); 



2. METHODOLOGY 

2.1. Estimator for non-Gausianity 

A measure of non-Gaussian properties inherent in a collec- 
tion of random samples is given by so called kurtosis statistic, 
(/C). For a set, S, of random samples, 5= {jt:,|/ = 1,2,3, ...,A^}, 
the /C statistic is defined as. 



J_ (x, -xo) 



(1) 



(119961): iBo uchet & Gispertl dl 99*91): 



Bennett et al. 

~T- 



legmark et al. (2003); Saha et al. (2006); Hansen et al. (2006); 
Eriksen et al. (2006); Hinshaw et al. (2007); Saha et al. ( 2008); 
Eriksen et al. (2008b a); Leach et al. (2008); Delabrouil le et all 
d2009i) ; ,Samal et al.. (.2010) : .Pietrobon et al.. (,20ia) ). For 
discussions about searching for non-Ga ussianity in fore - 
ground mini mized CMB maps we refer to Rath et al. (2009); 



where xo denotes sample mean and a is the standard deviation 
of samples. /C = for samples drawn from Gaussian distribu- 
tion whereas as discussed in Section 13.11 foreground distribu- 
tions possess large positive kurtosis. 

2.2. Foreground Cleaned Map 

Let X f denotes foreground contaminated CMB map in ther- 
modynamic temperature unit at 1 ° resolution at a frequency in- 
dex / and nb represents total number of available frequency 
bands. X/ is a column vector of n numbers where n is the num- 
ber of surviving pixels at a frequency band after masking the 
known point source positions. We form a foreground cleaned 
map following, 

IV, 

X'' = E^/X/. (2) 
/=i 

where Wf denotes the weight factor for linear combination for 
frequency index /. We choose these factors such that /C value 
for X"^ is minimized with an imposed constraint on them that 
Y^'f^iWf = 1. This condition preserves the actual values of 
CMB temperature at each pixel in the foreground cleaned map. 

Using Eqns. [T]and|2]one can show that, /C'^(W), of the lin- 
early combined map follows. 



Barreiro et al.,(,2000): Bandav et al.l(l2000): Raeth et al.i (,20ldl): 
Hou et alJ (l2010l) : lBernui & RebouqasI ( |2010|)~ 



/C"(W) = 



.^(WT,-W^ 



(3) 



WTW^) ,=1 

where W denotes a 1 x «/, vector whose /''' entry is given by 
Wf and T = Yl'j=\ T;- Here Tj is a «/, x «i symmetric matrix for 
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pixel j, Tjiff) = AffjATfij where Affj denotes the tempera- 
ture at pixel j of frequency index / after the mean temperature 
corresponding to this frequency has been subtracted from actual 
pixel temperature. 

In principle, the solution for W for minimum /C'^(W) satisfy- 
ing our conditions can be obtained by employing a Lagrange's 
undetermined multiplier approach. However, because of non- 
trivial nature of dependency of IC on W we find that such an 
approach is not feasible for our problem. Instead, we find the 
minimum of /C''(W) by invoking a non-linear search algorithm 
due to Powell. 

3. RESULTS 

3.1. Kurtosis 
I " Sy'nchro'tron,' k' banb' '__L_1' ' 
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Fig. 1. — Top Panel: Non-Gau.ssian distribution of pixel temperature of 
MEM foreground maps. Tlie frequency bands for this plots are chosen such 
that respective foreground component is most dominant over other frequencies 
inside the WMAP window. Bottom Panel: The plot of vaiiation of kurtosis, 
IC, as a function of foreground strength at WMAP frequencies. The sohd 
lines show mean curves obtained from the 150 random simulations of CMB 
and foreground with foreground strength as indicated by the ordinate of each 
points. The dashed lines show the variation of kurtosis for a randomly chosen 
realization. The error bars indicate cosmic variance. Mean kurtosis estimated 
from pure CMB maps is given by, (ICcmb) = —0.023 it 0.060. 

The distributions of pixel temperatures due to diffuse galac- 
tic foregrounds are shown in top panel of Fig. [1] Each of these 
foreground distributions is strongly asymmetric and exhibits 
a long tail towards positive temperature direction indicating a 
variation slower than ~ , which is the case for a Gaussian 
distribution. The peak for the synchrotron distribution shows 

^ Note that the position of known point sources are excluded before computing the 



that at 23GHz most likely contamination due to synchrotron oc- 
curs at a temperature ^ ISjiK, although contaminations at high 
(> AQQfiK) and low (e.g., as low as < ISjiK) pixel temperature 
are also likely. 

How does /C vary with foreground contamination in CMB 
maps at different WMAP frequencies? To answer this question 
we generate 150 random full sky CMB maps using WMAP's 
LCDM power spectrum. With each of these random realization 
of CMB map we add WMAP's MEM foreground templates at 
various strengths. The resulting behavior of the IC with fore- 
ground strength is shown in Fig. [T|for all WMAP frequency 
bands. For each band JC increases as the amount of foreground 
contamination increases. For K band the variation takes a shape 
of a plateau, characterized by a slow increase of JC, as the fore- 
ground level reaches ^ 20% of its full level. For all bands K, 
becomes more than 150 when all foregrounds are operative at 
their 100% level. It is interesting to note that the Q band (not 
the K band which has the highest level of synchrotron and free 
free contamination) shows the largest K, value among all the 5 
bands at the maximum foreground strength. 

3.2. Validation of the Method 

We validate the methodology by performing Monte-Carlo 
simulations. We use MEM foreground maps for synchrotron, 
free free and thermal dust available from LAMBDA website. 
These maps are provided in a common resolution of 1 ° and at 
pixel resolution parameter nside = 256 in antenna millikelvin 
temperature unit. We upgrade the pixel resolution of each map 
to nside = 512 and convert them to thermodynamic microkelvin 
unit. We add the composite foreground map of each WMAP 
frequency with a random realization of CMB compatible to 
LCDM model to make foreground contaminated CMB maps 
at each of the WMAP frequencies. Finally, we mask the posi- 
tion of the known point sources using the WMAP's 7 year point 
source mask. 

We develop a C code (hereafter GaussMap) to implement 
constrained Powell's conjugate gradient method. We perform 
200 Monte-Carlo simulations of foreground removal method. 
From these simulations we find that the foreground removal is 
effective over almost all parts of the sky. However, we find visi- 
ble signature of some residual foreground emission in the inner 
plane of the galaxy. To find out the sky regions where the resid- 
ual foreground could be significant compared to the expected 
CMB signal we subtract average of input CMB maps from the 
average of foreground cleaned maps. Then we form an initial 
mask by assigning zero values to all pixels with absolute tem- 
perature values more than 20/iK of this map and unity at all 
other pixels. This mask contains a set of scattered pixels in the 
inner galactic plane. To remove these pixels we first smooth 
the initial mask by a Gaussian window of 1°. We transform 
this smoothed mask to a new mask by assigning all pixels of 
smoothed mask with values greater than or equal to 0.9 to a new 
value of unity and the all other pixels to zero. To exclude posi- 
tion of the known point sources from the analysis we make the 
final mask multiplying this mask by the WMAP's point source 
mask. We call the resulting mask as G20 mask. This retains 
^ 87% of the entire sky area. 

We estimate JC for each of the cleaned maps obtained from 
Monte-Carlo simulations of foreground removal procedure Q. 
The mean kurtosis obtained from all the cleaned maps is given 
by {JC) = 0.86 ±0.21. This corresponds to ^ 4a detection of 

JC values. 



Removing foreground using non-Gaussianity 



3 



non Gaussianity from the cleaned maps. We interpret the non- 
vanishing (/C) in terms of residual foreground contamination 
originating from the galactic plane. After flagging of the pixels 
determined by the G20 mask we obtain, (K.) = -0.02 ± 0.06, 
which is consistent with zero. Henceforth, we use G20 mask 
as the basic mask to remove pixels contaminated by the resid- 
ual foreground while analyzing cleaned maps from Monte- 
Carlo simulations as well as WMAP data. The mean weights 
for 5 frequency bands satisfy, (W) = (0.049 ± 0.021 ,-0.419 ± 
0.063,-0.213 ±0.027, 1.643 ±0.063, -0.059 ± 0.029). 

3.2.1. Temperature Distribution 

Using the Monte-Carlo simulations we verify that the pixel 
temperature of cleaned maps outside the G20 mask follows a 
Gaussian distribution. For this we apply G20 mask both to a 
randomly chosen CMB realization and its foreground cleaned 
counterpart. To verify the Gaussian nature of these distri- 
butions we fit the histogram of input CMB map (after G20 
mask is applied) by a normalized Gaussian probability distri- 
bution, g{T) = exp(-{T- a)^ /(2s^))/ y^ilns^), where a and s 
denote respectively mean and standard deviation of the dis- 
tribution. From the fit we find that, s = 69.95 ± QW^iK and 
a = 2.12±0.1lAi/:B 

3.2.2. Power Spectrum 

We apply G20 mask on each of the 200 foreground cleaned 
CMB maps and estimate full sk y estimate of CMB p ower spec- 
ti-um using MASTER method (Hiv on et al.1 (l2002l) ). The av- 
erage of 200 power spectra obtained from the foreground re- 
moved maps matches excellently with the average of the input 
CMB power spectra. We show both spectra in Fig.|2]along with 
the cosmic variance. This verifies that no significant foreground 
contamination exists outside the G20 mask and the method out- 
lined in this paper can be used to estimate CMB power spectrum 
to extract cosmological information. 
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Fig. 2. — Comparison of average power spectrum obtained from fore- 
ground cleaned CMB maps outside G20 mask (red) witli the average input 
CMB power spectrum, (red line with eiTorbars). The average cleaned power 
spectrum is obtained by using MASTER algorithm. Both spectra have beam 
and pixel window effects removed. 



In some of the ea rlier pubUc ations (ISaha et al.l d200q) : 
ISouradeep etaLl(l2006l) : ISaha et al.l(l2008l) : ISamal et allOOlOl) ) 
the authors reported that the harmonic space based ILC method 

^ The small monopole is a manifestation of cosmic variance. 



gives rise to a negative bias in its power spectrum at the low 
multipoles. The bias appears due to a mere chance correlation 
between CMB and foregrounds at large angles due to availabil- 
ity of only small number of modes in the large scales on the 
sky. In principle, the negative bias in the power spectrum can 
be present in both pixel based and multipole based ILC algo- 
rithms which rely upon minimization of net foreground vari- 
ance from the cleaned map, since in these cases the finite cor- 
relation between CMB and foregrounds arises due to certain 
degree of overlap between hot spots and cold spots of CMB 
and foregrounds and the algorithm creates a negative bias due 
to nonli near dependence of weights on the empirical covariance 
matrix (ISaha et al.l (|2008)). However, the foreground removal 
method described in this paper relies upon the Gaussian na- 
ture of the final distribution, without explicitly minimizing the 
variance of the data. This leads to the advantage that the power 
spectrum obtained from the cleaned map does not have any neg- 
ative bias. 

3.3. Application on WMAP data 

Before the analysis we mask out the positions of known point 
sources from each of WMAP's 5 frequency maps using the 
point source mask. After masking each map contains 3054273 
pixels comprising 97% of the full sky area. 

3.3.1. Non-Gaussianity minimized map 
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Fig. 3. — Top panel: The foreground removed GMAP obtained by mini- 
mizing A'^(W). The map has a beam resolution of 1° and a pixel resolution 
corresponding to pixel resolution parameter nside = 512. Bottom panel: The 
difference between GMAP and the WMAP's 7 year's ILC map outside G20 
mask. 

We use GaussMap code to estimate the best fit 
weight factors corresponding to different frequency 
bands. The code gives best choice of weights as W = 
(0.035,-0.378,-0.217,1.632,-0.073) for K to W bands cor- 
responding to /Cm,„(W) = 1.17. The complete procedure to 
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estimate the weights takes only about a couple of minutes on 
an Intel 2.26 GHz processor. As one might expect, V band 
gets the maximum positive weight since it is supposed to be the 
least foreground contaminated frequency band in the WMAP 
observation window. The kurtosis for cleaned map outside the 
G20 mask becomes 0.093. We note that weights obtained from 
WMAP data have similar values to the mean weights obtained 
from Monte-Carlo simulations in Section IJ!2] 

Using the weights described above we obtain the foreground 
cleaned CMB map (top panel of Fig. [3] hereafter GMAP). Al- 
though there are some visible signature of presence of resid- 
ual foreground contamination in this map on the galactic plane 
(e.g. Cygnus region, Cassiopeia A, Carina nebula) all contami- 
nations are confined only near the galactic plane. We show the 
difference between GMAP and WMAP's 7 year ILC map out- 
side the G20 mask in bottom panel of Fig. [3] In the unmasked 
difference map larger pixel amplitudes are confined near the 
galactic plane, making a narrow strip-like structure, with ab- 
solute pixel temperature exceeding SOfiK. The bottom panel 
of Fig. [3] shows that applying G20 mask significantly reduces 
pixel amplitude of the difference map.. 



3.3.2. Temperature Distribution 
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Fig. 4. — The distribution of tlie pixel temperatures of GMAP outside G20 
mask is shown in red. The black curve shows a Gaussian fit to this data. The 
cyan curve shows the distribution of ILC map outside G20 mask. The other 
curves show the non-Gaussian distributions of pixel temperatures outside G20 
mask for each of the individual WMAP frequency band. Note that since these 
frequency maps contain CMB their distributions extend to negative pixel val- 
ues. 

We test for the distribution of pixel temperature of GMAP 
outside the G20 mask. The average pixel temperature of this 
map outside the G20 mask is only 0.85/i/r. Hence we fit for 
a Gaussian probability function with mean and variance s to 
the pixel temperature distribution of the masked GMAP. From 
the fit we obtain, s = 67.49 ±0.1 2/i/r. We show the distribution 
outside G20 mask obtained from GMAP and the WMAP's ILC 
map in Fig. H] We also show in this figure pixel temperature 
distributions of individual frequency bands. The long tails for 
K and Ka bands are due to strong synchrotron contamination. 
We note that all the distributions corresponding to 5 frequency 
bands are asymmetric representing their non-Gaussian nature. 
Pixel temperature of simulated frequency maps have distribu- 
tions similar to these histograms. 

3.3.3. Power Spectrum 



To estimate power spectrum from GMAP we first apply G20 
mask and estimate the partial sky CMB power spectrum. We 
convert the partial sky power spectrum to the full sky estimate 
by inverting the mode-mode coupling matrix. Finally we re- 
move both beam and pixel effect from this power spectrum. We 
show the resulting power spectrum in Fig. |5] As shown in this 
figure this power spectrum matches well with the ILC power 
spectrum estimated from the same sky region until / ^ 100. 
Beyond this multipole range GMAP power spectrum has less 
power than the ILC power spectrum. From the GMAP we find 
C2 = 246 AfiK^ and C3 = 4Q2.3fiK^ consistent with ILC esti- 
mates (C2 = 248. S/i/:^ and C3 = 404.0^*^:^). 
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Fig. 5. — Full sky estimate of CMB Power spectrum from GMAP outside 
the G20 mask is shown in red line. The full sky estimate of WMAP's ILC 
power spectrum computed from the same region is shown in blue hne. Both 
spectra have been connected for 1 ° beam effect and pixel window effect. In the 
inset we zoom in low multipole region. Both spectra match excellently with 
one another until / ~ 100. 

4. DISCUSSIONS & CONCLUSION 

In this paper we have developed, validated and applied on 
7 year WMAP data a global foreground minimization method 
from CMB sky based upon the measure of non-Gaussian nature 
of the diffuse galactic foregrounds. The GMAP obtained by 
this method matches well with the WMAP's ILC map outside 
the G20 mask. The power spectrum from GMAP also matches 
excellently with the power spectrum of ILC map until I ^ 100. 
At even higher / we find less power than WMAP's ILC map. 
The kurtosis value of the foreground minimized map outside 
the G20 mask, 0.093 ± 0.06, the error is estimated from Monte- 
Carlo simulations of foreground removal method. This is only 
a 1.5(7 effect implying that the CMB sky outside the G20 re- 
gion is sufficiently clean so that the signals from this region 
may be interpreted to have cosmological origin consistent with 
standard cosmological scenario. 

A crucial advantage of the method over the usual variance 
based ILC method is that the former does not possesses any 
negative bias at low multipoles. The quadrupole moment es- 
timated by us matches excellently with the WMAP's ILC es- 
timate. Our result shows that the problem of low quadrupole 
moment persist in WMAP data. The quadrupole and octopole 
maps also show similar nature as the ones estimated from the 
WMAP's ILC map. The problem of quadrupole and octopole 
alignment using GMAP would be investigated in a future pa- 
per - however, given the similarity between GMAP and ILC 
map (at the least, outside the G20 mask) the alignment likely 
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to remain. Another direction for future research would be to 
implement the method in the harmonic domain of the maps. 
The method can be generalized to estimate CMB angular power 
spectrum down to low angu lar sc ale followi ng approaches de- 
scribed in lSaha et al.l (l2006l) and ISaha et all 12008). Since po- 
larized foreground models are poorly known, an excellent re- 
search direction would be to apply our method on polarization 
sensitive CMB data released by WMAP and PLANCK satellite 
mission flTaub er (200i)). 

Our method shows that it is possible to remove foregrounds 
from the CMB sky purely based upon non-Gaussian nature of 
non-cosmological signals. However, if there is any cosmolog- 
ical signal of non-Gaussian origin the method would interpret 
them as foregrounds and hence would try to eliminate the sig- 
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nal from the cleaned map. Therefore comparing the GMAP 
with other foreground minimized maps which do not explic- 
itly rely upon the assumption of non-Gaussianity one may be 
able to isolate any possible primordial and/or secondary non- 
Gaussian signal. Estimating a Gaussian CMB map using data 
from current and future generation sensitive CMB experiments 
like PLANCK and CMBPol would play an important role to 
unfold this scientific information. 

We thankfully acknowledge useful discussions with Tarun 
Souradeep and Jayanta Kumar Bhattyacharjee. Some of the 
results of this pap er are o btained by using HEALPix soft- 
ware package ( Gorski et al.l l2b05) Fl We acknowledge use of 
LAMBDA data repository for our analysis. 
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